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Abstract 

We investigate by means of numerical simulations the possibilities of tomo- 
graphic techniques applied to a Bose-Einstein condensate in order to recon- 
struct its ground state. Essentially, two scenarios are considered for which 
the density matrix elements can be retrieved from atom counting probabili- 
ties. The methods presented here allow to distinguish among various possible 
quantum states. 
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I. INTRODUCTION 



Before the birth of quantum mechanics, optics and mechanics have long developed on 
parallel tracks, as light and massive entities were considered as waves and particles, respec- 
tively. At the beginning of the 20th century, with the introduction of quantum mechanics, 
waves and particles started to play an interchangeable role, with the concepts of photons and 
of De Broglie wavelength. This gave rise to the birth of quantum optics and atom optics [|I| . 
However, while an optical single-mode system has already been available since long ago, the 
same cannot be said for matter waves. In fact, in the field of atom optics, only recently 
breakthroughs in the evaporative cooling of dilute alkali gases have allowed the generation of 
Bose-Einstein condensates (BEC) 0. The BEC is a macroscopic occupation of the ground 
state of the gas and is one important paradigm of quantum statistical mechanics. 

In recent theoretical and experimental investigations of BEC, one of the most important 
and urgent issues is the determination of the actual quantum state of the condensate. In 
fact, as in optics, the presence of many particles in a single mode makes it possible to inquire 
about the multiparticle quantum state of the mode. At first thought, a number state might 
seem a natural description of a condensate mode, but the actual state may well depend 
on the details of preparation. For example the demonstration of first-order interference 
and observations of normalized, spatial correlation functions near unity suggest coherent 
states, while the presence of collisions between atoms may lead to the formation of squeezed 
states 1^. Also the internal states of the condensate atoms allow precise manipulation of 
the BEC state by interaction with light 

Few experimental methods of obtaining partial information about the state have been 
suggested 0. On the other hand, motivated by the success of quantum tomographic tech- 
niques in optics P), more direct methods for measuring the quantum state of a BEC were 
recently proposed p|-[TT|. 

In optical tomography, the key point is the use of a reference field, namely the local 
oscillator The latter, prepared in a coherent state, allows one to probe the desired state 
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through the measurement of a set of probabihties [0 . However, a difficulty arises if we try 
to adapt the same technique to a BEC. In fact, while in optics it is easy to obtain a coherent 
reference field (e.g. from a laser), the same is not actually available for atoms. Nevertheless, 
recent progresses in this direction seem promising Hence, in the present paper, we 

provide a detailed study of the possibility to reconstruct the quantum state of a BEC which 
includes both the scenarios: when a reference field is available and when it is not. 

The paper is organized as follows: in Section II we review the tomographic principle and 
consider a suitable operator transform on the atomic system. In Section HI, we consider 
the case of state reconstruction in the absence of a reference field, while in Section IV the 
opposite case is analyzed. Finally, in Section V we comment the numerical results, and in 
Section VI we conclude with a brief discussion. 

II. THE BASIC MODEL 



Recently, we established |]T5[ a quite general principle of constructing measurable prob- 
abilities, which determine completely the quantum state in the tomographic approach. For 
a more refined treatment see Ref. . 



Let us consider a quantum state described by the density operator p, which is a nonneg- 
ative Hermitian operator, i.e., 

P^ = P, Tr/5 = 1, (1) 

and 

{v\p\v)= p,,, > 0. (2) 

We label the vector basis | v) in the space of pure quantum states by the index v which 
may represent any degrees of freedom of the system under consideration. Formula (H) can 
be rewritten by using the Hermitian projection operator 

fl,=\v){v\, (3) 
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in the following form 

= Tr (n,p) . (4) 

On the other hand, in the space of states, there will be a family of unitary transformation 
operators U (cr) depending on some parameters a = (ai, . . . , 0"^ . . .), that can be sometimes 
identified with a group-representation operators. It was shown [|l5l that known tomography 
schemes can be considered from the viewpoint of group theory by using appropriate groups. 
More recently this concept has been developed obtaining an elegant group theoretical ap- 
proach to quantum state measurement ||16|. Here, we formulate the tomographic approach 
in the following way. Let us introduce a "transformed density operator" 

p. = U-\a)pU{a). (5) 

Its diagonal elements are still nonnegative probabilities 

{v \ Pa- \ v) = w{v, O") > . (6) 

These probabilities are functions of stochastic variable(s) v and parameter(s) a. As a conse- 
quence of the unit trace of the density operator, they also fulfill the normalization condition 

J dvw (f , a) = 1. (7) 

Of course, in the case of discrete indices, the integral in Eq. should be replaced by a sum 
over discrete variables. 

The l.h.s. of Eq. can be interpreted as the probability density for the measurement 
of the observable V (the operator whose eigenstates are given by 1^;)) in an ensemble of 
transformed reference frames labeled by the index a, if the state p is given. Along with this 
interpretation, one can also consider the transformed projector 

n,(or) = f/((T)n,[/-i(a), (8) 

in terms of which the expression (P) for the probability w {v, a) takes the form 

4 



w{v,a) =Ti[pIl,{a)] . (9) 
These probability densities are also called "marginal" distributions as a generalization of the 



concept introduced by Wigner |T7|. The tomography schemes are based on the possibility 
to find the inverse of Eq. (|^). If it is possible to solve Eq. considering the probability 
w {v, cr) as a known function and the density matrix as an unknown operator, the quantum 
state can be reconstructed in terms of measurable positive definite probability distributions. 
This is the essence of state reconstruction techniques. 

Specifically, we consider two atomic sources whose atoms (described by two bosonic 
modes bi and 62) can be mixed through an atomic beam splitter [Q, and assume that 
successively a phase shift (p can also be introduced between them. We shall specify these 
modes later. At the output, a detection of the number of atoms in both modes can be 
performed. This amounts to the possibility of measuring the probability distributions related 
to the transformed state 

p^U{9,^)pU^{9,4>), (10) 

where the transformation operator is given by 

U{e,(j)) =exp|-i^ [blb2exp{-i(f)) + biblexp{i(f))]^ . (11) 

Here, cos^(^/2) represents the transmission coefficient at the beam splitter. Eq. ( p!0D plays 
the same role of Eq. (^, and, in the spirit of the tomographic principle, the set of "marginals" 
associated to the transformed state will allow us to recover the original state. In the next 
two Sections, as anticipated in the Introduction, we shall distinguish two situations. 

III. CASE I 

We first treat the case where a reference field is not available. All we can do in this case 
is to consider two condensates belonging to the two modes 61, 62, and put the constraint of 
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total particle number conservation, i.e. [p, N] = 0, in order to infer their (joint) state. The 
latter is assumed to be a generic two-mode state of the type 

N 

\^) = J2cn\N-n),\n)2. (12) 

n=0 

At this stage we use the formal equivalence between the algebra for two harmonic oscillators 
and that for angular momentum |TP|. We write the state 11^)1 \N — 771)2 = \j + \j — 
as a spin state |m), where j = N/2 and m = n — j (m = —j, + — The j + 1 

states |m) have all the properties of the eigenstates of and Jz where 

j^ = jl=b\b2, l = l{b{h~bik) , = Jl + ]^{jJ. + J.J+) . (13) 

The effect of the beam splitter, including the phase shift, is a rotation by an angle 9 about 
an axis = u^- cos0 — Uy sin (p of the angular momenta J, i.e. 



U{e, 0) = exp (-i^J ■ u^j . (14) 
On the other hand, the rotation can be specified by means of the Wigner-D function 



{7n'\U{e, 0)|m) ^ p2^(V' = 0, 6, 0) , (15) 

where now if), 6, represent the Euler's angles. Then, the probability of j + m counts at 
the first detector and j — m at the second one is given by 

w{7n,9,c^)= j: f: 'D±^ii^,9,<P)pg\^^V^£:,^{7p,9,4>). (16) 

mi=-j m2=-j 

The measurement of the atomic number in both modes guarantees a unit efficiency. In 
fact, data for which the sum of counts is not can be disregarded. Moreover, in Eq. (|I6|) 
we have left the argument ip unspecified in the r.h.s. and omitted it in the l.h.s. since 
^mm' exp(— zm?/'): the marginal distribution only depends on the two angles 6 and 0. 
Following we will derive the expression for the density matrix of a spin state in 



terms of measurable probability distributions. This can be done by using the known integral 



product of three Wigner-D functions over the rotation group and the orthogonahty of the 
Wigner-3j symbols yV^^^^ms EH- Finally, the density matrix elements can be expressed in 
terms of the marginal distribution as 

2i j' 



j'=0 m,'=—j' m=—j 



X 0) WUL, W;^C,„, ^ (17) 

where the integration is performed over the rotation parameters, i.e. 

dn = J dtp J sin ed9 J dcj) . (18) 

Thus, Eq. ( plTl) can be used to sample two-mode BEC density matrix elements starting from 
the measurable probabilities w{m, 6, (p) and some known functions. 



IV. CASE II 



Recent progress in the generation of an atomic coherent source makes us hope about 
the possibility to have an atomic reference field M]. Thus, we shall consider in this section 
the first mode as the condensate to be investigated, and the second one coming from a 
coherent atomic source. 

Let p be the state of mode 1 we want to reconstruct and (3 the coherent state charac- 
terizing mode 2. Then, the probability of counting n atoms in mode 1 for = 7r/2, will 
be 



w{n,(3) = Tr U'^O = 7r/2, 0) p 



/22\ 



IJ{9 = 7r/2,0) \n)u{n\ 



(19) 



where (3 = \f3\exp{i(p), ip = aigjS — + 7r/2. This corresponds to the probability distri- 
bution for the measurement of the displaced number operator t)'^ {f3)h\hiD{(3) analogously 
to the Photon Number Tomography |^. In that case, however, one has to collect number 
distributions by spanning the whole complex plane /3; here, instead, we will simplify the 



procedure (see also @]). 



Of course, the number of atoms in the condensate, though not fixed, will be finite, thus, 
it happens that {k\p\m) = for k^m> Ni, with A^^i a suitable estimation of the maximum 
number of the condensed atoms. By virtue of this assumption we can rewrite Eq. (p!9| ) as 

w{n,P) = eM-\P\')nl E (fc|p|m)^==|/3r+'=-2" exp[z(m - fc)^]^-") 

xL(f-") , (20) 



where L^"^^ are the associated Laguerre polynomials. 

Let us now consider, for a given value of |/?|, the function w{n,P) as function of ip and 
calculate the coefficients of the Fourier expansion 

nj('\n,\(3\) = — f^" d^w{n,(3)exp{is^) , (21) 
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(s = 0, 1, 2, . . .). By combining Eqs. ( pOD and (|21|), we get 



Ni-s 

w^'\n,m = E Ai:lm){m + s\p\m), (22) 



m=0 



where 



J [171 + s)!m! 



If the distribution w{n, (3) is measured for n = 0,1,..., {N > Ni), then Eq. 
represents for each value of s a system of (A^ + 1) linear equations between the (A^ + 1) 
measured quantities and the (A^^i + 1 — s) unknown density matrix elements. Therefore, to 
obtain the latter we only need to invert the system 



N 



(m + s\p\m) = Y: M';:,lmW'\n, |/3|) , (24) 



n=0 



where the matrices A4 are given hj A4 = (A^A) ^A^. It is possible to see that such 
matrices satisfy the relation 

N 

E-^S,n(l/5|Mil(l/5|)= (25) 
n=0 

for m, m' = 0, 1, . . . , — s, which means that from the exact probabilities satisfying Eq. 
(^) the correct density matrix is obtained. By combining Eqs. (p^ and (0) we find that 
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(m + s\p\m 



ZTT nJ 

71=0 



(26) 



which may be regarded as the formula for the direct samphng of the condensate density 
matrix. In particular we see that the determination of the state of the condensate only 
requires the value of if (i.e. the phase between reference and condensate field) to be varied. 
Moreover, the presented reconstruction procedure involves Laguerre polynomials in place of 
additional summations,which guarantee a better stability in the numerical manipulation of 
large set of data, with respect to analogous methods . 

Finally, the non unit efficiency r] in the detection process can be accounted for by con- 
sidering a binomial convolution of the ideal probability p5| 



n=k 



(27) 



and the consequent modification of the matrix A. 



V. NUMERICAL RESULTS 

It is plausible, and it has been already suggested [p|,p!T|, that the state of a condensate 
with repulsive collisions be a squeezed state with reduced number fiuctuations. Hence in the 
following we will consider this situation. The single 1-mode state can be written as 

oo 

l^) = Ecn|n), (28) 

n=0 

where the coefficients Cn are given by 



where r is the squeezing parameter, xq is the (real) displacement and Hn denotes the Hermite 
polynomials. 



A phase space representation of this state can be given by the Q-function [|27| 



Q{a) = {a\p\a) , (30) 



where a is the complex amphtude of a coherent state. It yields 



Q{a) = exp(— |q;| 



oo 



E 



(31) 



n=o vn\ 

For the case discussed in Sec. Ill, we have to consider PJTT[] a two-mode squeezed state 
written in the angular momentum representation, i. e. 

j 



I*) E Cj+m\m) 



(32) 



m=-j 



where A/" is a normalization factor and the coefficients Cj+m are given in Eq. (pQ]). Since 
the quantity Xq + (r^ — l)/4r represents the mean number of atoms in mode 1, it must be 
smaller then the total number of atoms A^. Furthermore, the atomic coherent state basis 
for a system of angular momentum j is defined by ||2^ 



\9, 



j 

E 

m=-j 



( 



171 + j 



]-m 



sm - 

. 2, 



cos - 
. 2, 



exp(— zm0)|m) , 



(33) 



then, one can define the Q-quasiprobability distribution analogously to Eq. (^) 



Q(^,, 



; Ir n 5 



(34) 



For the state considered in Eq. 



it becomes 



m=-j 



(35) 



This is shown in Fig. 1(a). Instead, Fig. 1(b) displays the Q-function calculated from the 
reconstructed density matrix elements. We may see that the method of Sec. Ill is quite 
accurate, apart from some background noise. 

Analogously, in Fig. 2(a) we have plotted the ideal Wigner function of Eq. (pT]), while 
Fig. 2(b) is its reconstructed version. In this case, the statistical error depends on the chosen 
value of \(3\. For \f3\ close to zero the diagonal density matrix elements can be determined 
very precisely, whereas the off-diagonal elements strongly fluctuate. The opposite happens 
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by increasing \f3\. To compensate for the fluctuations, the number of measurement events 
must be increased. Another source of error stems from the truncation of the reconstructed 
density matrix at the value Ni. 

It is worth comparing the previous results with those obtained in the case of a number 
state. For this purpose, we show in Fig. 3(a) the ideal Wigner function for a number state. 
In this case, the Wigner function becomes negative, displaying the highly non classical 
character of a Fock state. As above, in Fig. 3(b) we show the Wigner function obtained 
by applying the reconstruction method of Sec. IV to such a state: the two pictures are 
practically indistinguishable, showing the accuracy of the present method. 

Finally, as an instructive comparison we show in Fig. 4 the same Wigner function of 
Fig. 2(a) calculated when the state reconstruction takes place with a random phase relation 
between probe and condensate. As can be seen the state becomes randomized and diffused 
in phase, but its Wigner function remains positive. This figure should be contrasted with 
Fig. 3. In this case the apparent U{1) symmetry is not pertaining to the state ||29|, but 
rather due to the measurement method, which implies a preparation of the (same) state at 
each experimental run. 



VI. CONCLUSION 

To conclude, we have studied, through numerical simulations, the possibilities of a tomo- 
graphic approach to the quantum state of a Bose-Einstein condensate. We have considered 
two possible scenarios, whether an atomic reference field is available or not. The correspond- 
ing methods turn out to be accurate and robust to detection inefficiency, and allow one to 
distinguish among various possible quantum states of the condensate. 

It is worth noting that the studied techniques allow direct sampling of the density matrix 



elements avoiding any ambiguities in the reconstruction procedure due to singularities 

The key point remains the possibility to have a reference field and/or its state prepara- 
tion. Furthermore, we note the necessity to deal with a relatively small number of atoms 
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in order to implement efficiently the numerical algorithms. In spite of these difficulties, we 
retain the possibility of measuring the true density matrix of a condensate accessible and 
worth considering. 

Finally, we would like to remark that the presented procedures could also be considered 
for other fields like high energy heavy ion collisions where pions can condense as well [ PT[] . 
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FIGURES 

FIG. 1. A squeezed 2-mode state, (a) Ideal Q-function when the displacement parameter is 
xq = \/5 and the squeezing parameter is r = e. (b) The corresponding Q-function reconstructed 
through the method of Sec. III. To obtain this figure we have simulated experimental data by 
adding to each probability w a noise term with a Gaussian distribution, the latter having a width 
proportional to the ratio between the probability itself and the number of runs for given parameters. 

FIG. 2. A squeezed state for the single mode, (a) Ideal Wigner function when the displacement 
parameter is xq = a/3 and the squeezing parameter is r = e. (b) The corresponding Wigner function 
reconstructed through the method of Sec. IV. The reconstruction parameters are = 1.1, rj = 0.9, 
and 3 x 10^ simulated experimental data per each phase have been used (see text). 

FIG. 3. A number state for the single mode, (a) Ideal Wigner function for the Fock state 
\n) = |5). (b) The corresponding Wigner function reconstructed through the method of Sec. IV. 
The reconstruction parameters are = 0.3, rj = 0.9, and 3 x 10^ simulated experimental data per 
each phase have been used (see text). 

FIG. 4. A squeezed state with a random phase between and 2it. Here, a value of the dis- 
placement parameter xq = a/5 has been used. 
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